Pre-dewetting transition on a hydrophobic wall: Statics and dynamics 
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For one-component fluids, we predict a pre-dewetting phase transition between a thin and thick 
low-density layer in liquid on a wall repelling the fluid. This is the case of a hydrophobic wall 
for water. A pre-dewetting line starts from the coexistence curve and ends at a surface critical 
point in the phase diagram. We calculate this line numerically using the van der Waals model and 
analytically using the free energy expansion up to the quartic order. We also examine the pre- 
dewetting dynamics of a layer created on a hydrophobic spot on a heterogeneous wall. It is from a 
thin to thick layer during decompression and from a thick to thin layer during compression. Upon 
the transition, a liquid region above the film is cooled for decompression and heated for compression 
due to latent heat convection and a small pressure pulse is emitted from the film into the liquid. 

PACS numbers: 64.70.F-,68.08.Bc,68.03.Fg 



I. INTRODUCTION 

Extensive efforts have been made on the wetting transi- 
tions for various liquids and walls both theoretically and 
experimentally |l|-|3|. As is well-known, when a liquid 
droplet is placed on a wall in gas, the three-phase con- 
tact angle changes from a finite value (partial wetting) to 
zero (complete wetting) at a wetting transition tempera- 
ture T w on the coexistence curve. Furthermore, there is 
a phase transition of adsorption between a thin and thick 
liquid layer across a prewetting line outside the coexis- 
tence curve [3, H-HI ■ For one-component fluids, the line 
starts from the wetting transition point T — T w on the 
coexistence curve and ends at a surface critical point at 
T = T ( F W in the T-n or T-p plane, where n and p are the 
density and the pressure in the surrounding gas region. 

On the other hand, many authors have been interested 
in the structural change in the hydrogen bonding net- 
work formed by the water molecules in the vicinity of a 
hydrophobic surface We also mention many ob- 

servations of surface bubbles or films on mesoscopic scales 
(with 10-100 nm thickness) on hydrophobic walls in wa- 
ter [l6-19]. Here the attractive interaction among the 
water molecules arise from the hydrogen bonding. As a 
result, there can even be a gas region in contact with a 
hydrophobic surface in liquid water at room temperature 
and at the atmospheric pressure (~ latom). To support 
this behavior, the solvation free energy AG so i of a hy- 
drophobic particle with radius R in water is nearly given 
by AttR 2 o per particle for relatively large i?(> lnm) [ll|, 
where a is the gas-liquid surface tension. 

In this paper, we examine whether or not a low-density 
film in liquid on a hydrophobic wall undergoes a pre- 
dewetting phase transion between a thin and thick low- 
density layer. As in the prewetting case, a pre-dewetting 
line should start from a point, T = Tf™, on the coex- 
istence curve and ends at a pre-dewetting critical point. 
It follows complete dewetting for T > Tj?™ on the coex- 
istence curve. We develop a mean-field theory based on 
the Ginzburg-Landau model as in the original paper [l[ . 
Our calculations are thus performed rather close to the 



critical point (at T ~ 0.9T C ), where the film density is 
not very small compared to the ambient liquid density. 

Furthermore, we are interested in the dynamics of 
the wetting transition of volatile liquids, where under- 
standing of evaporation and condensation at the inter- 
face is still inadequate [2fj| - r22j . For example, Koplik et al. 
[23j performed molecular dynamic simulation to observe 
evaporation of a droplet and a decrease of the contact 
angle upon heating a wall in partial wetting. Guena et 
al. [24j performed an experiment, where a weakly volatile 
droplet spread as an involatile droplet in an initial stage 
but disappeared after a long time due to evaporation in 
complete wetting. In a near-critical one-component fluid, 
Hegseth et al. [25] observed that a bubble was attracted 
to a heated wall even when it was completely wetted by 
liquid in equilibrium. To study such problems, we have 
recently developed a phase-field model for compressible 
fluids with inhomogeneous t emp erature, called the dy- 
namic van der Waals model [26j | . (S ee a review on vari- 
ous phase-field theories of fluids j27|.) In our framework, 
we may describe the gas-liquid transition and convective 
latent heat transport without assuming any evaporation 
formula. We then numerically investigated evaporation 
of a liquid droplet on a heated substrate [28[ and spread- 
ing of a liquid film on a cooled or warmed substrate for a 
one-component fluid [29| . The lattice Boltzmann method 
has been applied to two-phase fluids also (30l - [32l |. How- 
ever, this method has not yet been fully developed to 
describe evaporation and condensation. 

This paper also presents simulation results on the film 
dynamics u sing the numerical method in our previous 
studies [ID, l29j | We initially start with an equilibrium 
thin or thick film on a wall at the bottom near the pre- 
dewetting transition and then cool or heat the tempera- 
ture at the top. (We use "bottom" and "top" though we 
do not assume gravity.) Subsequently, the cell is gradu- 
ally decompressed or compressed and the pre-dewetting 
transition is induced in the film. 

The organization of this paper is as follows. In Sec. II, 
we will examine the static aspect of the pre-dewetting 
transition in the Ginzburg-Landau scheme [l], 0] in the 



mean field theory. In Sec. Ill, we will apply the dynamic 
van der Waals model to investigate the pre-dewetting dy- 
namics by cooling and heating the top plate of a cylin- 
drical cell. In the appendix, the pre-dewetting transition 
will be examined near the critical point by expanding 
the free energy with respect to the density n around the 
critical density n c up to the quartic order. 



II. STATICS 

A. Ginzburg-Landau model 

We consider a one-component fluid in contact with a 
sold wall in equilibrium, where the number density n is 
the order parameter. Assuming short-ranged forces, we 
set up the free energy F with the gradient contribution 
as HE SI 



F= dr 



f(n,T) + ±M\Vn\ 



+ / da-yn, (2.1) 



where the integral is in the fluid container in the first 
term and on the wall surface in the second term ( J da 
being the surface integral). As a function of n and T, 
f = f(n,T) is the Helmholtz free energy density. In 
our numerical_calculation, we will use the simple van der 
Waals form |H, 



/ = k B Tn 



ln[nA t 3 h /(l - v n)} - 1 



CVqTI 



(2.2) 



where Vq is the molecular volume, e is the magnitude of 
the attractive pair potential, and A t h = h(2n /mksT) 1 ' 2 
is the thermal de Broglie length with m being the molec- 
ular mass. The coefficient M of the gradient free energy 
will be assumed to be independent of n but proportional 
to T. The last term is the surface free energy expressed 
as the integration on the solid surface [lH3|, where the 
second-order term of the form A , /A~ 1 (n — n c ) 2 /2 (present 
in the original work [l|) is neglected. (This is allowable 
for very large |A e | compared to the correlation length 
[35l].) In the literature [3 GUI, the so-called surface field is 
given by —7. For water-like fluids, 7 > for a hydropho- 
bic surface and 7 < for a hydrophilic surface. For this 
surface free energy, the pre-dewetting transition occurs 
for 7 > 0, while the prewetting transition for 7 < 0. 

In equilibrium, the space-dependent density n — n(r) 
in the bulk region is determined by 



/' 



MV 2 n 



Mo 



(2.3) 



where p = df/dn is the chemical potential and fig is 
a constant. The left hand side fj, — MW 2 n represents 
the generalized chemical potential including the gradient 
contribution. On the wall surface we have 



Mv ■ Vn = 7, 



(2.4) 



where v is the outward normal unit vector (from the wall 
to the fluid) on the wall surface. 



B. Pre-dewetting transition for no > n c 

Let us consider an equilibrium one-component fluid in 
the region z > in contact with a planar substrate with 
7 > placed at z = 0. In this one-dimensional geom- 
etry, all the quantities depend only on z. The density 
n(z) tends to a constant liquid density no far from the 
wall. The pressure and the chemical potential far from 
the wall are written as po and /io, respectively, where 
Pa = nok B T/(l — vqUq) — voeriQ for the van der Waals 
model (2.2). In equilibrium, we should minimize the 
grand potential f2 (per unit area) given by 



O = 



dz[u 



-M\n' 
2 1 



7"s 



where n' = dn/dz and n s is the surface density, 
n s = n(0). 

We introduce the grand potential density as 

w = / - non + po, 
Use of the van der Waals form (2.2) gives 

7l(l - VqUq) 



(2.5) 



(2.6) 



(2.7) 



= k B Tn\n 



n (l - v n) 
ev (n - n Q ) 2 . 



-k B T- 



n 



1 - v a n 



Far from the wall, w behaves as u> 
and tends to zero. We define 



(2.8) 

x(no)(n - n ) 2 / 2 



x (n) = d 2 f/dn 2 = l/n 2 K T , 



(2.9) 

n(z) is 



where Kt is the isothermal compressibility. 

In the region z > 0, the density profile n 
determined by \i— Mn" — (io from Eq.(2.3), where n" = 
d 2 n/dz 2 . We multiply this equation by n' and integrate 
the resultant one with respect to z to obtain [33[ 



lu = M \n'\ 2 /2. (2.10) 

From Eq.(2.4) the boundary condition at z = reads 

Mn'(0)=7. (2.11) 

On the basis of the van der Waals model (2.2), we intro- 
duce a microscopic length t defined by 

£ = (M/2k B Tv ) 1/2 . (2.12) 

Since 7 > 0, Eq.(2.10) is rewritten as 

n'(z) = v / 2w(n(z))/M = l~ 1 y/u}(n(z))/k B Tv , (2.13) 

where lo is treated as a function of n and its T dependence 
is suppressed. Using £, we integrate this equation as 



= £ 



dn/^/uj(n)/k B Tv . 



(2.14) 



CO 



CM 




surrounding fluid approaches a state on the coexistence 
curve. In our case, the film thickness D increases log- 
arithmically as no approaches the liquid density on the 
coexistence curve n^ x , where we may take the isothermal 
path of approach, for example. To show this, we expand 
the free energy density / around the corresponding gas 
density on the coexistence curve n^ x as 



/ 



-Pc 



/'c 



Xg(t 



(2.17) 



where p cx and u cx are the pressure and the chemical 
potential, respectively, on the coexistence curve. When 
the bulk density no is slightly larger than n^ x , we have 
bo - Pc X )/n e cx = fi - u cx = xe( n o - "-ex) i so that 



FIG. 1: (Color online) ^/2Muj(n)/e£ vs v n at T/T c = 0.9 
for vono = 0.5525, 0.56,0.58,0.60 from below, where the 
smallest one is very close to the coexistence liquid density 
n cx — 0.5524. Use is made of the van der Waals form (2.8). 



w = Xg (n - + Xl (n^ - < x )(n - r&). (2.18) 

Substitution of this relation into Eq.(2.14) yields 

n s ) 



D = L\n 



^Xg( n cx 



Xi(n 



(2.19) 



From Eqs.(2.11) and (2.13) we have 

7 = ^2Mu(n s ) = 2£ v / k B Tv uj{n s ). (2.15) 
Using Eq.(2.13) we rewrite Eq.(2.5) as [l[ 



J are defined on the 
I 1 / 2 is the correlation 
tends to zero, a 



where the integrand vanishes at the lower bound n = n s 
from Eq.(2.15). 



0.08 



0.05 



0.02 



where x g = x(«? x ) and Xe = x( n L> 
coexistence curve and £ g — (M/Xc X ; 
length in the gas phase. As n — n l c 

well-defined interface appears. 

In Fig.l, we plot y/2Moj(n)/et = 2y/k B Tv uj(n) /e vs 
von at T = 0.9T C for four bulk densities no, using the van 

der Waals model (2.8). (i) The smallest bulk density is 

£1=1 dn[^2Mu(n) - 7] + -fn , (2.16) 0.5525u ~ 1 , which is slightly larger than the coexistence 

liquid density n^ x = 0.5524i' ~ 1 . As a result, uj(n) nearly 
vanishes at the coexistence gas density ng x = 0.1419u,7 1 . 
The gas layer thickness D is logarithmically dependent on 
n — n e cx as in Eq.(2.19). (ii) At n — 0.56W0 -1 , there are a 
minimum and a maximum satisfying duj/dn = f — uo = 
0. (iii) At no = 0.5854WQ -1 , the two extrema merge into 
a point n = n sc , at which d 2 w/dn 2 = X — holds and 
n sc coincides with the so-called spinodal density on the 
gas branch for the linear form of the surface free energy, 
(iv) At no = 0.60VQ 1 , ^Ju{n) increases with decreasing 
n in the range n < no- Then Eq.(2.15) yields a unique 
surface density n s for any 7 > 0. 

As illustrated in Fig. 2, a pre-dewetting transition ap- 
pears when the curve of y/uAn) has two extrema as in 
the case of n — 0.56VQ 1 in Fig.l. That is, if the areas 
Si and S2 are equal, a first-order transition occurs be- 
tween two surface densities n^ and . For each given 
7 > 0, a pre-dewetting line starts from a temperature, 
T = 1^(7), on the coexistence curve and ends at a pre- 
dewetting critical temperature, T — T^ w (^), outside the 
coexistence curve. The wall is completely dewetted by 
the gas phase for T > T^(-f) on the coexistence curve. 

In Fig. 3, two examples of the pre-dewetting line are 
written for j/ei = 0.05 and 0.1 in the T-n plane (left) 
and in the T-p plane (right), where n and p denote those 
in the bulk (no and po)- In this paper, the coefficient M 
is independent of n and is proportional to T as 
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F IG. 2: I llustration of a pre-dewetting transition, where 
^/2Mui{n)/el is plotted for T/T c = 0.860 and v n = 0.59407. 
At 7/e^ = 0.05, the areas of the regions Si and 52 coincide, 
leading to the same value of Q in Eq.(2.16) at two surface 
densities n s = n, and n? . 



Cahn [l[ showed that the thickness of a film of the 
preferred phase on a wall grows logarithmically as the 



M = TC. 



(2.20) 
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FIG. 3: (Color online) Pre-dewetting lines (pdwl) for j/el = 
0.05 and 0.10 in the T-n plane (left) and in the T-p plane 
(right), where n and p are the bulk values (written as no 
and po in the text). Each pdwl starts from a point on the 
coexistence curve, where T = T^(-y) and p — p^ij), and 
ends at a pre-dewetting critical point, where T = T dw (7) 
and p = Pc(l)- See Sec. Ill for nonequilibrium simulations 
starting with points A and B as initial equilibrium states. 
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FIG. 4: (Color online) Left: (a) Pre-dewetting temperature 
T^c W on the coexistence curve and (b) pre-dewetting critical 
temperature T dw as functions of j/el. Their linear approx- 
imations are (a') T c d x w '° and (b') T dw '° (see Eqs.(2.22) and 
(2.23)). Right: (a) Pre-dewetting pressure p d ™ on the coex- 
istence curve and (b) pre-dewetting critical pressure p dw as 
functions of "i/el. 



Then C and hence I = (C/2wofcs) 1 / 2 are indepen- 
dent of n and T. However, essentially the same re- 
sults were obtained even if M is independent of T 
(not shown in this paper). We give T cx w and T dw 
and the corresponding densities njr™ and n dw . For 
7 /e£ = 0.05, we have (T c d x w /T c , v n^) = (0.851,0.601) 
and (T dw /Tc,v < w ) = (0.943,0.524). For 7 /e£ = 
0.1, we have (T™/T c , v n^) = (0.683,0.724) and 



(T dw /T c ,v n dw ) = (0.883,0.606). 

In Fig. 4, the left (right) panel displays the pre- 
dewetting transition temperature (pressure) T^(j) 
(Pcxil)) on the coexistence curve and the pre-dewetting 
critical temperature (pressure) T dw (7) (p dw (7)) as func- 
tions of j/el. For each 7, the pre-dewetting line is in 
the range T d x w ( 7 ) < T < T dw ( 7 ). If j/el < 1, these 



> 




0.85 0.87 0.89 0.91 0.93 0.95 

T/T c 
Surface density 



o 
> 



vono=0.59 




y/6i=0.05 



T/T c 



FIG. 5: (Color online) Normalized desorption wo^ _1 rid in 
Eq.(2.25) (top) and normalized surface density von s (bottom) 
as functions of T/T c for vono = 0.59, 0.57, 0.55, and 0.53 (from 
left) at "//el = 0.05. They change discontinuously across the 
pre-dewetting transition. Theoretical curve for vo£~ 1 Tid from 
the Landau expansion is shown for vono = 0.53 (broken line, 
rightest in top plate). 



temperatures are both close to T c and are expanded with 
respect to 7 as 



T dw ( 7 ) 



T c (l 
T c (l 



A cx7 /e£ +•••), 
A cl /el +••■)• 



(2.21) 
(2.22) 



On the basis of the van der Waals model (2.2), the coef- 
ficients A cx and A c are calculated in the appendix as 



9 



2^+3, A c 



(2.23) 



We also plot the linear approximations T c d { W:0 (7) = T c (l — 
A c ^/e£) and T dw '°( 7 ) = T c (l - A c j/e£), which are in- 
deed in good agreement with the numerical curves for 
small 7 <C el. For the usual prewetting transition in 
the case 7 < 0, the counterparts of T c dw and T dw are the 
wetting temperature T w on the coexistence curve and the 
prewetting critical temperature T c pw . As will be discussed 
in the appendix, they are expanded as in Eqs. (2.21) and 
(2.22) if 7 is replaced by |7|. 

In Fig. 5, we show the desorption T\d per unit area in 
the top plate and the surface density n s in the bottom 
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FIG. 6: (Color online) Normalized density profile von(z) vs 
normalized distance z/l from a wall at y/e£ = 0.05. Shown are 
those for T/T c = 0.93665, 0.93664, and 0.92 at v n = 0.52892 
(top) and those for v n = 0.54358, 0.54357, and 0.52892 at 
T/T c — 0.92 (bottom). The first two curves from above are 
taken just before and after the pre-dewetting transition, while 
the third one close to the coexistence curve. 



C. Rough estimates of C and 7 for water 

In our continuum theory, the constant C in Eq.(2.21) 
and the length £ in Eq.(2.12) remain arbitrary. For each 
fluid, we may roughly estimate their sizes with input 
of the experimental values of T c , n c , and the surface 
tension a at some temperature T. For example, we 
have T c = 647K and n c = 1.076 x 10 22 cm~ 3 for wa- 
ter. Using these T c and n c in the van der Waals model 
calculation, we have e = 28k B T c /8 = 3.0 x 10" 17 mJ 
and «o = (3n c ) _1 = 31 x 10~ 24 cm 3 , from which we 

1 /3 

define the van der Waals radius by a v( jw = v — 
3.1 A. For water, an experimental surface tension is 
(7 = 9.3mJ/m 2 = 0.030e/a 2 at T — 0.92T C [3^. On 
the other hand, from the model (2.1) under Eqs.(2.12) 
and (2.21), the surface tension is numerically calculated 
as er v dw = 0.0097e£/«o at T = 0.92T C . If we equate these 
experimental and numerical values of the surface tension, 
we obtain £ = 3.1a v dw and C/2k B ~ 10 a vdw S3]- ^ n tne 
same manner, for argon, we have e = 0.70 x 10 _17 mJ, 
v = 4.11 x 10 21 cm 3 , and 7 = 9.3mJ/m 2 = 0.016^ at 
T = 0.92T C H, leading to C/2k B - 6.80^. 
For water, we estimate the surface field 7 by 

7 - (r/nuq, (2.25) 

for a typical hydrophobic surface. Here nn q and a are 
some appropriate liquid density and surface tension. We 
use the the above-mentioned surface tension for water 
and set nn q ~ n c to obtain 7/e^ ~ 0.03. This estima- 
tion is based on the behavior of the solvation free energy 
AG so i of a large hydrophobic particle in water [l5| . as 
discussed in Sec.l. 



plate as functions of T/T c for vouq = 0.59, 0.57, 0.55, and 
0.53 at 7 1 el — 0.05. These lines are outside the coex- 
istence curve in the T-n plane. They change discontin- 
uously across the pre-dewetting transition in the range 
T c d x w < T < T c dw . The discontinuities vanish as T — > T c dw 
We calculate the desorption of the fluid by 

/>oo 

Ti d = / dz[n Q -n(r,z)}. (2.24) 
Jo 

On each line, as T T c d x w , T ld grows logarithmically as 
(no — n 9 cx )D with D being given by Eq.(2.20), while n s 
tends a well-defined limit (see the lower panel of Fig. 5). 
The curve for voUq — 0.53 is closest to the pre-dewetting 
criticality, so we also display its theoretical approxima- 
tion from the Landau expansion in the appendix. 

In Fig. 6, we display the density profile n(z) at "f/e£ = 
0.05. In the top plate we set T/T c = 0.93665,0.93664, 
and 0.92 at fixed uq ~ 0.52892v ~ 1 , while in the bottom 
plate we set v n = 0.54358, 0.54357, and 0.52892 at fixed 
T = 0.92T C . In these plates, the first two curves repre- 
sent the profiles just before and after the pre-dewetting 
transition, while the third one is obtained close to the 
coexistence curve with a well-defined thick gas layer. 



D. Nanobubbles on a hydrophobic spot 

Dewetting as well as wetting is very sensitive to het- 
erogeneity of the substrate. Here we realize bubbles in 
equilibrium on a heterogeneous surface with a position- 
dependent 7. 

In this paper, we suppose a circular hydrophobic spot 
with radius 25^ on the bottom surface by setting 

7(r) = 0.05e£ (r < 25f), 

= (r > 25£), (2.26) 

where r = (x 2 + y 2 ) 1 / 2 at z = 0. The resultant equi- 
librium density n(r, z) satisfies Eq.(2.3). Its boundary 
condition at z = is given by dn/dz = O.Ob'j/M = 
0.025e/k B Tv £ for r < 25£ and dn/dz = for r > 251 

Slightly above the pre-dewetting line, the top plate of 
Fig. 7 gives the cross-sectional density profile n(r, z) at 
y = for the point A in Fig.3, where T = 0.94T C , p = 
0.89p c , and n = 0.52894WQ" 1 . The minimum density is 
n s = 0.40587f; _1 at the spot center on the surface. Below 
the pre-dewetting line, the bottom plate of Fig. 7 gives 
the profile at the point B in Fig.3, where T — 0.92T C , 



Equilibrium outside pre-dewetting line 




Equilibrium inside pre-dewetting line 



and heat and mass fluxes passing through the interface 
remain weak. As a result, T and the generalized chemical 
potential (the left hand side of Eq.(2.4)) are continuous 
across the interface. 

In our diffuse interface method, the interface thickness 
needs to be longer than the simulation mesh size Ace, so 
our system size cannot be very large. In our simulation to 
follow, its length is 150^. Nevertheless, our cell contains 
many particles about Vno ~ 6 x 10 6 ^ 3 /urj, where V is the 
volume and no is the initial liquid density. 




FIG. 7: Equilibrium axisymmetric density profile n(z,r) at 
y — on a hydrophobic spot with radius 25£, within which 
j/e£ = 0.05 and outside which 7 = 0. Top: T = 0.94T C , p = 



0.. 



and no = 0.52894t> at point A in Fig. 3, where the 



fluid is outside the prewetting line and the minimum surface 
density is n s — 0.40587i; ~ 1 on the center of the spot. Bottom: 
T = 0.92T C , po = 0.71177p c , and n = 0.52925U0 -1 at point B 
in Fig.3, where the fluid is inside the prewetting line and the 
minimum surface density is n B — 0.088433HO -1 . 



Po = 0.71177p c , and no = 0. 52925^ . The minimum 
density is n s = O.O88433V0 -1 at the spot center. These 
profiles are very different. To characterize the film size, 
we introduce the desorption of the fluid by 



2tt 



dz 



drr[no — n(r, z)], 



(2.27) 



where the gas film is well within the integration region 
< z < zq and r < tq. By setting zq = 15£ and r$ — 
37.5^, we obtain V = 339^0 1 and 3460£ 3 w " 1 for the top 
and bottom plates in Fig. 7, respectively. 



III. DYNAMICS 

In this section, using the dynamic van der Waals 
modelprjj, we numerically investigate the dynamics of 
a thin gas layer created on the hydrophobic spot in 
Eq.(2.27). We treat a one-component fluid without grav- 
ity in a temperature range 0.92 < T/T c < 0.94, where 
the gas density is 25 — 30% of the liquid density (see 
the bottom curves in Fig. 6). Then the mean free path 
in the gas £ m { is not very long. For very long £ m f, how- 
ever, numerical analysis based on a continuum phase-field 
model becomes very difficult. In our simulation, further- 
more, relatively small temperature changes are applied 



A. Hydrodynamic equations with gradient stress 

We set up the hydrodynamic equations for the mass 
density p — mn, the momentum density pv, and the 
entropy density S [H,[2§|, where m is the molecular mass 
and v is the velocity field. Here S consists of the usual 
entropy density ns and the negative gradient entropy as 



S = ns-ic|Vn| 2 . 



(3.1) 



where s is the entropy per particle and the coefficient 
C in Eq.(2.21) appears here. The internal energy den- 
sity, written as e, can also contain the grad ient contri- 
bution, but we neglect it for simplicity [2g|. Then the 
total Helmholtz free energy density is given by e — TS = 
/ + M|Vn| 2 /2 as in the first term as given in Eq.(2.1). In 
our scheme, we use the entropy equation instead of the 
energy equation to achieve the numerical stability in the 
interface region. That is, with our entropy method, we 
may remove the so-called parasitic flow at the interface, 
which has been encountered by many authors (ioj . 

We integrated the following hydro dynamic equations 
without gravity for p, pv, and S [28l[29|: 



+ V • (pv) = 0, 

^-pv + V • (pvv + S) = V • a, 
at 



01 



S + V 



Sv - Cn(V • v)Vn 
+ {e v + e e )/T, 



V • -VT 

T 



(3.2) 
(3.3) 

(3.4) 



where the terms in the right hand sides are dissipative. 
In Eq.(3.3), II = {Ily} is the reversible stress tensor, 



n, 



p- CT(nV 2 n+ -|Vn| 2 ) 
-CT(V i n)(V i n), 



6 U 



(3.5) 



where p is the van der Waals pressure. Hereafter Vj = 
d/dxi with Xi representing x, y, or z. The terms propor- 
tional to C arise from the gradient entropy, constituting 
the gradient stress tensor. The a — {<Jij} in the right 
hand side of Eq.(3.3) is the viscous stress tensor, 



ri(y iVj 



(3.6) 



in terms of the shear viscosity rj and the bulk viscosity £. 
In Eq.(3.4), A is the thermal conductivity, while 



e v = yVijVjtfr) e'e = A(VT) 2 /T, 



(3.7) 



are the nonnegative entropy production rates arising 
from the viscosities and the thermal conductivity, respec- 
tively. On approaching equilibrium, e v and eg tend to 
zero, leading to vanishing of the gradients of v and T. 
The (total) energy density in the bulk is defined by 



eT 



pv 2 /2, 



(3.8) 



which includes the kinetic energy density. The energy- 
conservation equation reads [39| 



d 



-V- 



e T v + (I! - a) ■ v - AVT 



(3.9) 



If Eqs.(3.2) and (3.3) are assumed, the entropy equation 
(3.4) and the energy equation (3.9) are obviously equiva- 
lent. In the text book [39(, however, the entropy equation 
is derived from the three fundamental conservation equa- 
tions (3.2), (3.3), and (3.9) Now the total fluid entropy 
iStot an d the total fluid energy £ tot are written as 



Stot = / drS, £ to t = / dre T + J dwyn. (3.10) 

If the surface field 7 is independent of T as in Eq.(2.26), 
the surface term in Eq.(2.1) is the surface energy and 
there is no surface entropy on all the boundaries. In a 
fixed cell, we assume the no-slip condition v = on its 
boundaries. Then the space integrations of Eq.(3.4) and 
(3.9) yield the time derivatives of <Stot and £tot- 



d f e v + e e 



d fl "' AW + 7 ", (3.11) 



T 



d f 

—Sua = / da{v AVT + 7 n), 



(3.12) 



where h — dn/dt and 7 can be heterogeneous as in 
Eq.(2.26). If there is no heat input and n = on the 
boundary walls, <5>tot continue to increase monotonically 
until an equilibrium state is realized. 



B. Simulation method 

We suppose a cylindrical cell in the axisymmetric ge- 
ometry, where our model fluid is in the region < z < H 
and < r = (x 2 + y 2 ) 1 ^ 2 < L. Assuming that all the 
variables depend only on z, r and t, we perform simula- 
tions on a two-dimensional 300 x 300 lattice. We take the 
simulation mesh length Ax equal to Ax =1/2, where I 
is defined in Eq.(2.12). Thus our cell is characterized by 



The velocity v vanishes on all the boundaries. The vis- 
cosities and the thermal conductivities are proportional 
to n as 



r) = C = vomn, A = 4fceJ/on. 



(3.14) 



These coefficients are larger in liquid than in gas by the 
density ratio ne/n g (~ 5 in our simulation). The kine- 
matic viscosity vq — rj/mn is a constant. We will mea- 
sure time in units of the viscous relaxation time, 



r = I 2 l /v = C/2k B v Q u , 



(3.15) 



H = 1504 L = 150£, 



(3.13) 



on the scale of t. The thermal diffusion constant £>th = 
X/Cp is of order vq, where C p is the isobaric specific 
heat per unit volume of order ksn (not very close to 
the critical point). The time mesh size At in integrat- 
ing Eqs.(3.2)-(3.4) is 0.02ro. If the dynamic equations 
are made dimensionless, there appears a dimensionless 
number given by a = rav^/el 2 (written as a in Ref.[26]), 
where m is the molecular mass. The transport coeffi- 
cients are proportional to vq oc ct 1 / 2 . In this paper we 
set a — 0.06, for which sound waves are well-defined as 
oscillatory modes for wavelengths longer than I (2(| . 

As the boundary conditions, we assume Eq.(2.4) for 
the density n (even in nonequilibrium) so that Cdn/dz = 
0.05e^/T on the hydrophobic spot at z — and v-Vn = 
on all the other surface regions. For the velocity v, the 
no-slip condition v = is assumed. The temperature 
is fixed at To at the bottom z = and at Th at the 
top z — H. The side wall is thermally insulating as 
dT/dr = at r = L. 



C. Decompression by cooling the upper plate 

We prepared the equilibrium state at T = 0.94T C in the 
upper plate of Fig. 7 as an initial state at t = 0. We then 
cooled the top temperature T H from 0.94T C to 0.92T C 
fixing the bottom temperature To at 0.94T C for t > 0. 
After a very long time (t > -ff 2 /4T> t h), the fluid tended 
to a steady heat-conducting state. 

In Fig. 8, we show the profiles of the temperature 
T(r, z,t) at r = and the density n(r,z,t) at r = L/2 
as functions of z to illustrate how cooling and decom- 
pression are realized in the cell. In the initial stage, the 
piston effect i s op erative, which has been studied the- 
oretically [Jl], |42[ and experimentally [43Tj45l |. In this 
situation, there appear a compressed thermal diffusion 
layer at the top and an expanded one at the bottom with 
thickness growing in time as ^d(^) = (Dtht) 1 ^ 2 , which 
produce sound waves propagating in the cell and causes 
an adiabatic change outside the diffusion layers. The 
sound velocity in the present case is about c = M/t , so 
the acoustic traversal time over the cell is H/c ~ 35ro- 
In fact, at t = 500to, T and n are flat in the middle re- 
gion, where T and n are changed adiabatically. However, 
in the late stage t > IOOOtq, the layer thickness reaches 
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FIG. 8: (Color online) Normalized temperature T(r, z,i)/T c 
at r = (top) and density von(r,z,t) at r/L = 0.5 (bot- 
tom) vs z/H after cooling the top from 0.94T C to 0.92T C . At 
t = 500t~o, thermal diffusion layers appear near the top and 
bottom, causing an adiabatic change in the middle region. For 
t > lOOOro, the thermal diffusion extends throughout the cell. 
The density is decreased in the lower part of the cell, result- 
ing in a homogeneous pressure decrease. Around t = 1320ro, 
a cool spot appears above the film during the pre-dewetting 
transition from a thin to thick film. In the inset, the profiles 
of T and n at r — and L/2 near the bottom are compared. 
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FIG. 9: Normalized temperature T(r, z,t)/T c (top) and pres- 
sure p(r, z,t)/p c (bottom) as functions of t in units of To in 
Eq.(3.15) at (z/H,r/L) = (0.1,0.1) (slightly above the gas 
film) after cooling the top. Two insets in each panel show ini- 
tial oscillatory behavior due to sound wave propagation (left) 
and a sudden change at the prewetting transition (right). 



H/2 and the thermal diffusion becomes relevant through- 
out the cell. In all these processes, the pressure p is kept 
nearly homogeneous in the cell. 

In our problem, we should focus on the behavior of 
T and p around the gas film on the hydrophobic spot. 
The film is under gradual decompression, but nearly at 
the initial temperature (since the bottom temperature is 
pinned). As a result, at t — 1320Yo, the pre-dewetting 
transition takes place from a thin to thick gas film. This 
is indicated by formation of a cool spot in the liquid 
above the film in Fig. 8. It is caused by the latent heat 
adsorption to the expanding film at the first-order phase 
transition (see Fig. 12 in more detail). There is no cooling 
outside the film region r/L > 0.2. 

In Fig. 9, we display the temperature T(r, z, t) and the 
pressure p{r, z,t) slightly above the film at a fixed posi- 
tion located at (z/H, r/L) = (0.1, 0.1). In the early stage 



t < 300to, we can see their oscillatory relaxations caused 
by sound wave traversals. Upon occurrence of the pre- 
dewetting transition at t ~ 1300to, we can see a small 
drop in T of order 5 x 10 _4 T C and a small peak in p of 
order 10 _3 p c on their curves. However, the insets in Fig. 9 
reveal that their behavior is somewhat complicated on a 
short time scale of order 50to, because they change on 
arrivals of a sound wave and a velocity disturbance from 
the film. At long times, T tends to a constant about 
0.938T C , but p continues to decrease slowly. 

In Fig. 10, the profiles of the density n(r,z,t) vs z are 
given at r/H = 0.1 at three times together with their 
cross-sectional profiles. The increase in the film thick- 
ness is abrupt around t — 1320to at the pre-dewetting 
transition and is very slow at t — 2000to. Figure 11 
presents the time evolution of the surface density n s (r, t) 
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FIG. 10: Normalized density von(r, z,t) vs z/H at r/ff = 0.1 
for t/ro = 0, 1320, and 2000 from above after cooling the 
top. At t = 1320to, the film thickness is increasing abruptly 
at the pre-dewetting transition, while at t — 2000i~o it is still 
increasing slowly (see Fig. 11). Cross-sectional film profiles 
are also given in the right, where the darkness represents the 
density at four times. 
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FIG. 11: Normalized excess desorption iio^ _3 Ar(t) in 
Eq.(3.16) and surface density von s (r,t) vs t at r/L = 0.1 
in units of to after cooling the top. They exhibit oscillatory 
behavior in the initial stage and an abrupt change at the pre- 
dewetting transition taking place in a time interval of 250ro . 



at r/L — 0.1 and the excess desorption Ar(t) denned by 

rZo pro 

Ar(i) = 27r/ dz drr[(n (r, z) -n(t,r,z))], (3.16) 
Jo Jo 

where no(r, z) = rt(0, r, z) is the initial density profile and 
we set zq — 0.1H and ro = 0.25L. The initial oscillatory 
relaxations arise from traversals of sound waves, while n s 
decreases and T increases abruptly at the pre-dewetting 



FIG. 12: (Color online) Temperature T(r,z,t) in gradation 
at t/ro = 1000, 1320, and 2000 around the hydrophobic spot. 
On the block line, |Vn| is largest, indicating the film location. 
Arrows indicate the velocity v. Their amplitudes are known 
from reference arrows below the middle abd bottom panels. 
The liquid region above the film is cooled due to latent heat 
absorption to the film at the pre-dewetting transition. The 
profile at t = 2000ro (bottom) is nearly stationary, where 
condensation is taking from the side and evaporation from 
the upper surface. The colors represent T/T c according to 
the color bar at the bottom. 



transition at t ~ 1.3 x 10 3 Tn. We recognize that the pre- 
dewetting transition from a thin to thick film occurs in a 
time of order 250r . 

In Fig. 12, we show T(r,z,t) in gradation at t/ro — 
1000, 1320, and 2000 around the hydrophobic spot. Since 
there is no clear interface here, a gas-liquid boundary 
is indicated by a line on which |Vn| is largest in the 
direction of Vn. We also display the velocity v by 
arrows. The reference arrow below each panel repre- 
sents the maximum velocity, being equal to 0.0363, 1.48, 
and 0.0878, at i/<r = 1000,1320, and 2000, respec- 
tively, in units of 10~ 2 i?/Tn. It is enhanced in the middle 
plate at t/ro — 1320 during the pre-dewetting transi- 
tion, where we can see an upward flow with a magni- 
tude of 1.5 x 10~ 2 ^/ro in the liquid region above the 
expanding film. The expanding velocity of the film is 
of the same order. If we multiply this velocity by the 
duration time 250to of the transition (which is inferred 
from Fig. 11), we obtain a film thickness of order D ~ At 
in accord with the density profiles in Fig. 10. The cor- 



responding Reynolds number around the film is about 
D\v\/vq ~ 0.1. Remarkably, the region above the film 
is cooled due to latent heat absorption to the film by 
2 x 10 _3 T C . The profile at t = 200to is nearly stationary, 
where the typical velocity is of order 5 x 10 _4 £/to. A 
balance is attained between condensation from the side 
and evaporation from the upper surface, while the fluid 
is at rest far from the film in the presence of a con- 
stant temperature gardient. Finally, we give the heat flux 
Qb(r, t) — —AdT/dz at z = from the wall to the fluid 
at r = and L/2: (Q b (0,t),Q b {L/2,t) = (1.39,1.13) at 
t = 1000, (5.71,0.942) at t = 1320, and (0.618,0.867) 
at t = 2000 in units of \Q- A ei/v T . The heat flux is 
much enhanced at the center during the pre-dewetting 
transition. In the inset of the upper panel of Fig. 8, the 
gradient \dT/dz\ is very large at r = 0, but it is much 
smaller at r = L/2. 



D. Compression by heating the upper plate 

We prepared the equilibrium state at T = 0.92T C in the 
lower plate of Fig. 7 (at the point (B) in Fig. 3) as an initial 
state at t = 0. We then heated the top temperature Tjj 
from 0.92T C to 0.94T C with the bottom temperature T 
held fixed at 0.92T C for t > 0. 

In Fig. 13, the time evolution of T(r, z, t) at r — and 
n(r 7 z 7 t) at r — L/2 is illustrated as functions of z. As 
in Fig. 8, the piston effect takes place in the initial stage, 
but in the reverse direction. That is, at t = 500to, there 
appear an expanded thermal diffusion layer at the top 
and a compressed one at the bottom. In the late stage t > 
1000to, the thermal diffusion extends over the cell. The 
pressure p is nearly homogeneous in the cell and increases 
slowly in time above the initial value. Thus the film 
is under gradual compression, but nearly at the initial 
temperature. Occurrence of the pre-dewetting transition 
is indicated by formation of a heat spot in the liquid 
above the film at t = 1280to. It is caused by the latent 
heat release from the shrinking film (see Fig. 17). 

Figure 14 displays the temperature T(r,z,t) and the 
pressure p(r, z, t) at (z/H,r/L) = (0.1,0.1) as in Fig. 9. 
The adiabatic process takes place in the early stage. We 
find occurrence of the pre-dewetting transition at t <~ 
1280to, where T and p exhibit a small peak of order 2 x 
10~ 4 T C and 2 x 10~ 3 p c , respectively, with a duration time 
about 50to. Their detailed behavior can be seen in the 
insets of Fig.14. At long times, T tends to 0.9216T C , but 
p continues to increase slowly. 

In Fig. 15, the profiles of n{r,z,t) vs z are given at 
r/H = 0.1 at three times together with their cross- 
sectional profiles. The film thickness decreases gradu- 
ally around t = 1000to and is stationary at t = 2000to. 
Figure 16 gives the time evolution of the surface density 
n s (r,t) at r/L = 0.1 and the excess desorption AT(t) 
defined in Eq.(3.16). Here, n s (r,t) increases abruptly 
around t = 1000to as in Fig.ll. However, the decrease in 
AT(t) is steep in the early stage t < 200tq and is gradual 
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FIG. 13: (Color online) T(r,z,t)/T c at r = (top) and 
von(r, z,t) at r/L = 0.5 (bottom) vs z/H after heating the 
top from 0.92T C to 0.94T C . At t = 500to, thermal diffusion 
layers appear near the top and bottom. For t > 1000ro, the 
thermal diffusion extends throughout the cell. The density is 
increased in the lower part of the cell, resulting in a homo- 
geneous pressure increase. Around t = 1280ro, a heat spot 
appears above the film during the pre-dewetting transition. 
In the inset, the profiles of T and n at r = and L/2 near 
the bottom are compared. 



later until t ~ 1.3 X 10 3 To. Evaporation from a thick film 
is significant in the early stage before the transition. 

In Fig. 17, we show T(r,z,t) in gradation at </tq = 
500, 1280, and 1500 around the film. The film location 
is indicated by a line on which |Vn| is largest along Vn. 
The velocity v is displayed by arrows. We can see a down- 
ward flow from the liquid region to the film induced by 
the film shrinkage. The reference arrow below each panel 
represents the maximum velocity, being equal to 0.284, 
3.54, and 0.0215 at t/r = 500, 1280, and 1500, respec- 
tively, in units of H)~ 2 1/tq. The shrinking velocity of the 
film is of the same order. In contrast to the cooling case 
in Fig. 12, the region above the film is somewhat heated 
even in the early stage (at t/ro = 500), which is because 
of the film shrinkage in Fig.ll. The heating is most en- 
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FIG. 14: T(r,z,t)/T c (top) and p(r, z, t)/p c (bottom) as func- 
tions of t in units of To at (z/H,r/L) — (0.1,0.1) (slightly 
above the film), after heating the top. Two insets in each 
panel show initial oscillatory behavior due to sound waves 
(left) and a sudden change at the pre-dewetting transition 
from a thick to thin film (right). 
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FIG. 15: v n(r,z,t) vs z/H at r/L = 0.1 for t/r = 0, 1000, 
and 2000 from below after heating the top. The film thickness 
decreases as in Fig. 16. It is stationary at t — 2000ro. Cross- 
sectional film profiles are also given in the right, where the 
darkness represents the density at four times. 
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hanced at t/ro — 1280 in the middle plate during the pre- 
dewetting transition. The profile at t = 1500rn is nearly 
stationary, where a balance is attained between evapora- 
tion from the side and condensation from the upper sur- 
face. Finally, we give the heat flux Qj,(r,t) = —XdT/dz 
at z = from the wall to the fluid at r = and L/2. That 
is, (-Q b (0,t),-Q b (L/2,t) = (1.24,1.22) at t = 500, 
(44.2,0.870) at t = 1280, and (1.02,0.871) at t = 1500 
in units of 10 _4 e£/i>nTn. The heat flux is much enhanced 
at the center during the pre-dewetting transition. In the 
inset of the upper panel of Fig. 13, we can see that dT / dz 
is large at r — and is much smaller at r = L/2. 



IV. SUMMARY AND REMARKS 

We summarize our main results, 
(i) In Sec. II, for one-component fluids, we have ex- 
amined the pre-dewetting transition on a hydrophobic 
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FIG. 16: (Color online) Normalized excess desorption 
VqI~ Ar(t) in Eq.(3.16) and surface density von s (r,t) vs t 
at r/L — 0.1 in units of to after heating the top. Here n s 
changes abruptly during the pre-dewetting transition in the 
time region 800ro < t < 1300ro, but AF(i) decreases steeply 
in the initial stage t < 200ro. The system becomes stationary 
for t > 1300r . 



wall in the mean-field theory. We have assumed the 
Ginzburg-Landau free energy for the number density n 
in the bulk and the surface free energy linear in n with 
7 representing the interaction between the wall and the 
fluid. Depending on the sign of 7, a line of prewetting 
or pre-dewetting appears in the gas or liquid side of the 
coexistence curve. If 7 is small, the line is close to the 
critical point as in Eqs.(2.21) and (2.22). In numerical 
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FIG. 17: (Color online) Temperature T(r,z,t) at t/ro = 
500, 1280, and 1500 around the hydrophobic spot. On the 
block line, |Vnj is largest. Arrows indicate the velocity v. 
Their amplitudes are known from reference arrows below the 
panels. Even the early stage, the film thickness decreases and 
the region above the film is heated due to latent heat release. 
The profile at t = 1500tq is nearly stationary. 



analysis in Figs. 3,5, and 6, we have set 7 = 0.05e£ 
(where e is denned in eq.(2.2) and £ in Eq.(2.12)). For 
water a rough estimation yields such an order of magni- 
tude of 7. For real systems it is desirable to measure 7. 
In Subsec.IID, we have numerically realized a localized 
film on a hydrophobic spot in Eq.(2.26), where a thin 
film is above a pre-dewetting line at T = 0.94T C and a 
thick one is below it at T = 0.92T C . 
(ii) In Sec. Ill, we have investigated the time evolution of 
a localized film on the hydrophobic spot in the axisym- 
metric geometry. Starting with the equilibrium films 
in Fig. 7, we have numerically integrated Eqs.(3.2)-(3.4) 
after a change of the top temperature from 0.94T C to 
0.92T C in Subsec.IIIC and from 0.92T C to 0.94T C in Sub- 
scc.IIID. Cooling at the top leads to a pressure decrease, 
while heating at the top leads to a pressure increase. 
We have found singular behavior at the pre-dewetting 
transition in the hydrodynamic variables and the excess 
desorption AT(i) in Eq.(3.16). Most markedly, the liquid 
above the film is cooled for decompression and is heated 
for compression due to latent heat convective transport 



from the growing or shrinking film, as shown in Figs. 8, 
12, 13, and 17. A small pressure pulse is emitted from 
the film to propagate through the cell, as shown in Figs .9 
and 14. This pulse could be detected experimentally (45T |. 

We give some remarks. 

(1) A pre-dewetting line as well as a prewetting line 
readily follows near the critical point for a small surface 
field 7. It is then of great interest where we can find a 
pre-dewetting line for water on a given hydrophobic wall 
in the phase diagram. 

(2) In future work, we should examine the role of the 
long-ranged van der Waals interaction 0-11] on the 
pre-dewetting phase transition. 

(3) For water-like polar fluids, a small amount of 
impurities can strongly promote phase separation due 
to the solvation effect [46], though we have treated 
one-component fluids only. As a solute, we may add 
a noncondensable gas (such as CO2), hydrophobic 
particles, or ions in water. Such impurities can strongly 
affect the formation of surface bubbles or films on a 
hydrophobic wall in water [16l4l9l ] . We will report shortly 
on how the pre-dewetting line is shifted downward with 
increasing the gas concentration. 

(4) If the mesh length Ax = ijl is a few A, our system 
length is on the order of several ten manometers and the 
particle number treated is of order 10 7 — 10 8 (see the 
beginning of Sec. III). Our continuum description should 
be imprecise on the angstrom scale. Thus examination 
of our results by very large-scale molecular dynamics 
simulations should be informative. We should also 
investigate how our numerical results can be used or 
modified for much larger film sizes. 

(5) Phase changes inevitably induce a velocity field 
carrying heat and mass. It has been crucial during the 
pre-dewetting transition. We have found a steady flow 
around the film at long times in Figs. 12 and 17. In our 
previous simulation [291 ] . a steady circular liquid film 
was realized on a homogeneous wall in the complete 
wetting condition under a critical heating rate, where 
evaporation and condensation balanced. It evaporated 
to vanish for stronger heating, while it expanded for 
weaker heating or for cooling. 

(6) We should study the two-phase hydrodynamics 
such as evaporation or boiling, where the hetrogeneity 
of the wall is crucial, as exemplified in this paper. In 
fluid mixtures, a Marangoni flow decisively governs the 
dynamics with increasing the droplet or bubble size even 
at very small solute concentrations (47j . 
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Appendix: Calculations near the critical point 

We examine the pre-dewetting transition near the crit- 
ical point in the mean-field theory. The surface free en- 
ergy density is linear in n as in Eq.(2.1). The density n 
and the temperature T are assumed to be close to their 
critical values. We expand the Helmholtz free energy 
density / = f(n,T) in powers of the density deviation 
%j) = n — n c up to the quartic order as 

f(n.T) = f c (T) + u c (T)ip + ^ 2 + (Al) 

where f c (T) = f(n c ,T), fj. c (T) = fi(n c ,T), and 

r = T/T c - 1 (A2) 

is the reduced temperature. We assume r < 0. Then 
the liquid and gas densities on the coexistence line are 
n c + ip e and n c — ip e , respectively, with 

^ e = ( ao |r|M>) 1/2 . (A3) 

In the van der Waals model in Eq.(2.2) [H, T c = 
8e/27fcs, n c = v^ 1 /3, and 

a = 2ev 0l u = 9ev%/2. (A4) 

With the Landau expansion (Al), the grand potential 
density u> in Eq.(2.7) is written as 

u=^W>- + + 3^ 2 - 2</> e 2 ), (A5) 

where ipQ — uq — n c is the value of ip far from the wall 
assumed to be small. The derivative w' = du;/dn is equal 
to the chemical potential deviation fi — fio written as 

J = u (^ - V'oXV' 2 + M + - €)■ (A6) 
Here we introduce a dimensionless surface field 7 by 

7 = 7(2/ Uo M) 1 / 2 ^7 2 

= 9Vay/8e£|T|, (A7) 

where the second line is the result of the van der Waals 
model. For each 7 and ipo/ip e , the condition (2.15) yields 
the equation for the surface density deviation ip s — n s — 
n c (which is smaller than ip ) in the form, 

(V>o - + + 3V 2 - 2^ 2 ) 1/2 /^e = 7- (A8) 

As illustrated in Fig. 2, this equation has two solutions 
tpg = rig — n c and ipf = — n c , corresponding to the 
two surface densities at the pre-dewetting transition. 

First, we take the limit no — > n% x or ipo —> ip e on 
the prewetting line, where the film thickness grows as 



in Eq.(2.20) but tp^, ip^ , and 7 tend to finite limiting 
values. From Eq.(A8) we obtain 7 = \ip 2 B /tp 2 — 1] in this 
limit. Since ips < — V'e an d ips > V'e we find 

^=n L s -n c = -V, e (l+7) 1/2 , 
i>» = n»-n c = ^ e (l-7) 1/2 . (A9) 

The equal area requirement S\ — S2 in Fig. 2 yields 

(l+ 7 ) 3 / 2 -(l- 7 ) 3 / 2 = 2, (A10) 

whose square gives 37 2 — 1 = (1 — 7 2 ) 3 / 2 . Again taking 
the square of this equation, we find 7 + 67 2 = 3, which 
is solved to give 

7= \J2VZ - 3 = 0.68125. (All) 

Therefore, as T — > 7^(7), the two surface densities at 
the pre-dewetting transition are ips — — 1.2966^/> e , and 
ipf = O.5646-0 e . From Eq.(A7), the pre-dewetting value 
of 7 on the coexistence curve is written as 

7c X (T) = ^2/V3-l) 1 / 2 |rM. (A12) 

We then obtain the coefficient A cx in Eq.(2.24). 

Second, we seek the value of 7 at the pre-dewetting 
critical point, denoted by j c (T). Since d 2 ui/dn 2 — at 
ri = n s , n s coincides with the spinodal density on the gas 
branch so that 

ip s =n s -n c = -tp e / V3. (A13) 

Then from Eq.(2.15) we find 

7c (T) = 8\r\e£/9. (A14) 

From the additional condition to' = at ip = ip s , Eq.(A7) 
yields the bulk density deviation ipo in this case as 

ipo = n - n c = 2^i e /v3, (A15) 

which is larger than ip e as it should be the case. For 
each small 7, the reduced temperature r and the den- 
sity deviation no — n c are given by r^ w = — 9^/8e£ and 
4| T c w |/3v / 3^0) respectively, at the pre-dewetting critical 
point. We then obtain the coefficient A c in Eq.(2.24). 

Our calculation results are applicable to the prewetting 
transition near the critical point, where 7, ipo, and tp s are 
small negative quantities. To use the above relations, we 
should set ijj e = — (ao|r|/uo) 1//2 and replace 7 by |7|. See 
the comment below Eq.(2.23) on T w and TP W . 

We should note that Papatzacos [48| obtained the wet- 
ting angle 9 W in partial wetting for the linear surface free 
energy in Eq.(2.1) and the bulk free energy in Eq.(Al) 
(see Ref.[30l] also). In our notation, the dewetting angle 
#dw = 7T — W satisfies 

(1 + 7 ) 3/2 - (1 - 7) 3/2 = 2 cos0 dw , (A16) 



where the fluid is on the coexistence curve. The above 
equation becomes 7 6 + 6j 4 + 3j 2 {1 — 2cos(2#d w )] = 
sin 2 (2#dw)- Setting (3 = 3 _1 arccos(sin 2 #dw), we solve 
this equation to obtain 

7 = 2sgn(7r/2 - 6 dw ) Vcos/3(l - cos/3), (A17) 



where sgn(x) gives the sign of x. For 7 > 0, we have 
the complete dewetting limit #dw — > (or 6 W — > n) as 
T -> T c d x w = T c (l - A^jel + ■ • •) from below. For 7 < 0, 
we have the complete wetting limit 9 W — > (or 6dw — > i") 
as T -> T w = T c (l - AacM/e^ + ■ • •) from below. 
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